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Frequency combs with weakly lasing exciton-polariton condensates 
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We predict the spontaneous modulated emission from a pair of exciton-polariton condensates due 
to coherent (Josephson) and dissipative coupling. We show that strong polariton-polariton inter¬ 
action generates complex dynamics in the weak-lasing domain way beyond Hopf bifurcations. As 
a result, the exciton-polariton condensates exhibit self-induced oscillations and emit an equidistant 
frequency comb light spectrum. A plethora of possible emission spectra with asymmetric peak dis¬ 
tributions appears due to spontaneously broken time-reversal symmetry. The lasing dynamics is 
affected by the shot noise arising from the influx of polaritons. That results in a complex inhomo¬ 
geneous line broadening. 

PACS numbers: 42.55.Ah, 78.67.-n, 71.36.+C, 42.79.Hp 


Condensation of exciton-polaritons (EP’s) in semicon¬ 
ductor microcavities formed by two distributed Bragg 
mirrors with quantum wells between them has been ex¬ 
perimentally observed hhsi- Being incoherently excited 
in the microcavity, EP condensates are in general out of 
thermodynamic equilibrium. EP condensates refuel their 
particle depot through absorption of cavity photons and 
emit coherent light due to tunneling of the composite 
EP states through distributed Bragg mirrors. Sample 
inhonrogeneity, either accidental or intentional, can in¬ 
duce several condensation centers (CC’s) (3j [6HS] ■ At low 
enough pumping, one expects a system of disconnected 
BEC droplets emitting light at different uncorrelated fre¬ 
quencies. As the pumping increases the condensates tend 
to establish mutual coherence and emit in a laser mode 
|6j. Already two CCs can synchronize and emit at a single 
joint frequency mug. This is possible because the con¬ 
densates exchange particles due to Josephson coupling 
and adjust their emission frequencies, which in turn de¬ 
pend on the number of condensed particles due to the 
polariton-polariton repulsion. In addition to the coher¬ 
ent Josephson coupling, there can be a dissipative (ra¬ 
diative) coupling between CC’s, which reflects the de¬ 
pendence of the losses in the system on the symmetry of 
singe-particle states. That new stationary regime called 
weak lasing takes place when pumping rates reside be¬ 
tween some minimal and maximal rates of losses [10] . In 
the weak lasing regime, the system is stabilized by the 
formation of specific many-particle states which adjust 
the balance between gain and loss in the system. 

In this Letter we show that in the weak lasing regime 
two CCs can emit not only at a single frequency, but 
also at a whole frequency comb which in principle con¬ 
tains an infinite number of equidistant lines of coherent 
lase-like radiation. This emission reflects the fact of for¬ 
mation of spontaneous selfsustained anharmonic oscilla¬ 


tions of both the occupation numbers and the relative 
phase between the condensates, in sharp contrast to pre¬ 
viously reported damped Josephson oscillations HUES- 
We study possible emission spectra and the way they are 
affected by noise. While the emission frequency of single- 
line EP lasers resides in the eV range mmmmm, the 
modulation frequency of comb emission can be adjusted 
to be in the terahertz and sub-terahertz range. Filtering 
out of the high-frequency component through optical de¬ 
modulation yields the low-frequency coherent signal as a 
new promising type of coherent terahertz emitters. The 
EP self-induced oscillation is also a novel mechanism of 
optical frequency comb generation as compared to mode- 
locked lasers mm and optical microresonators [HI HE 
Consider two coupled EP condensates with order pa¬ 
rameters 

Vh,2 = y/n&e**^, (1) 

where nj .2 are the occupations of the two condensates, $ 
is the total phase and 2</> is the phase difference. The time 
evolution of ^ 1,2 is governed by the Langevin equations 

(n= i) pm 

- - -A/v + aj-i/vl 2 ^) + / M (t), (2) 

where n ^ v = 1,2 label the condensates . The pa¬ 
rameter g = r — W describes the difference between the 
rates of losses T and pumping W, denote the singe- 
particle energies of the condensates, the parameters 7 
and J define dissipative and coherent coupling between 
the condensates, respectively, and a is the polariton- 
polariton interaction constant. The last term in Eq. © 
is the Gaussian white noise satisfying = 0 and 

(/mW/^'A')) Due to gauge invariance, 
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only the frequency detuning ui is relevant and in what fol¬ 
lows we will count the frequency from ujq = (oq + oj 2 )/2. 
Rescaling time we can fix 7 = 1 and, since rescaling the 
condensate amplitudes is equivalent to a change of a, we 
can set a = 2 without loss of generality. 

The dissipative coupling induces a relative phase ip de¬ 
pendent dissipation in the system. This can be observed 
from the eigenvalues A which control the condensate evo¬ 
lution ipi t 2 ~ e At in the absence of interaction [TO]. For 
sufficiently large 7 jg one of the eigenmodes turns unsta¬ 
ble. Therefore, the dissipative coupling acts as a phase- 
selective pump which depends on the relative phase <p\ it 
pumps one eigenmode while keeping the other one lossy. 
In this regime nontrivial weak lasing states are formed 
(see Ref. [10] for a complete account). 

First we consider the noise-free case (f\p = 0) in 
order to determine the attractors of the system (Hi- 
Two nontrivial solutions F ± to Eqs. ([ 2 ]) were identified 
in Ref. [1(1] , They are characterized by nonzero time- 
independent triplets {ni,n 2 ,ip}± with the total occupa¬ 
tion n = ni+712 = [(g 2 + J 2 )R ± +gui]/gJR ± , n\ = n(l — 
JR ± )/2, n 2 = n( 1 + JR ± )/2, <j) = (n — arctani? ± )/2, 
where R± = ±^/(l — g 2 )/(g 2 + J 2 ). The total phase $ 
of the condensates satisfies 


= 



J cos(2 (j>)\ 
4,/ni n 2 ) 


(ni + n 2 ) + 


sin(2</>) 

4 y /n 1 n 2 


(ni - n 2 ), 


(3) 


and the rhs in ([ 3 ]) gives a time-independent frequency 
= —STo for these solutions. The two centers evolve 
in a coherent fashion 7 / 7 2 ~ e - *^ 0 * and ft 0 dehnes the 
blue-shift of the emission line with respect to the average 
single-particle frequency. In the subspace {ni,n 2 ,<j)}± 
these states correspond to fixed points [19]. In parts of 
the control parameter space these states are stable, and 
should manifest themselves as weak lasing states. The 
F ± states loose stability at g = gf [T9 


2 (at + = (at + 1 ) [w 9c + (at + j 2 )R ± ] (4) 


We plot the two instability curves in the { g , w} space at 
fixed J = 0.1 in Fig. [T] (solid lines). The F ± states are 
unstable in the shaded areas LC 1 * 1 . In particular they are 
both unstable in the joint area LC + and LC - , where the 
trivial solution n\ = n 2 = 0 is unstable as well. What are 
then the stable stationary states of the system, if any? 

The answer is obtained by linearizing the phase space 
flow around F ± in the subspace {ni, n 2 , 4>}±- At g = 
g'lf two corresponding eigenvalues are purely imaginary 
±iAft, with their real parts changing sign. As a result, a 
supercritical Hopf bifurcation occurs, where stable limit 
cycles LCA with frequency 

Aft = y/2 g 2 + J 2 /g 2 + J 2 + ujg/R ± (5) 


are born around the respective unstable fixed points F 1 * 1 

HDH2H- 



shaded areas in the w, ^-parameter space. Stable LC 1 * 1 are 
born through Hopf bifurcations at the solid line and turn 
unstable at the dashed gf lines where they undergo period 
doubling bifurcations. They are the only stable attractors to 
coexist in the central (yellow) region. Here J = O.I 7 . 


Away from the bifurcation line the LCs increase the os¬ 
cillation amplitudes, deform, and change their frequency. 
The coexistence region of LC + and LC - grows in size as 
the Josephson tunneling is reduced. At the Hopf bifur¬ 
cation, where a LC emerges, ni j2 ,ip and also $ become 
periodic functions of time with period T = 2tt/ Aft. Then 
they may be expanded in a Fourier series with frequency 
harmonics iVAft and N = 0,±1,±2,... The integration 
of the constant term (N = 0) in the Fourier series of $ 
results in a linear time dependence, $dc = similar 

to F ± . Therefore 

VvW =P^)e~ inot i (6) 

where the functions p M (t) = p^,(t + T) are periodic in 
time. The Fourier spectrum of ^(t) is equidistant with 
frequency harmonics positioned at fto + A/Aft. 

Approaching the dashed lines gf in Fig0 the cor¬ 
responding LC turns unstable and undergoes a period 
doubling bifurcation. This gives rise to a new stable 
period-doubled LC, which however again quickly un¬ 
dergoes a period doubling bifurcation. A period dou¬ 
bling route to chaos along a Feigenbaum scenario leads 
to chaotic attractors [21]. Therefore just two coupled 
exciton-polariton condensates suffice to produce an ex¬ 
tremely rich and complex synchronized dynamics. 

Experimentally the polariton order parameter is de¬ 
tected by analyzing the emitted light from the microcav¬ 
ity. In near-field measurements, only small parts of the 
sample, like one condensation center, can be probed. Our 
aim is to calculate the spectral density 7^2 (ft) of the ra¬ 
diation corresponding to different nontrivial attractors. 
Applying a Fourier transformation (FT) we have 

4(ft) = |FT(^(t))| 2 , 


M = 1,2. 


(7) 
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FIG. 2: (Color online) Asymmetric frequency combs in the 
near-field spectrum from a LC + (a),(b) and a period doubled 
LC + (c). Here w = 0, J = 0.1. The small arrows indicate 
the position of the respective N=0 peak, (a) I\ with g = 0.19 
(upper panel), g = 0.27 (lower panel), (b) The same as (a) but 
for 1 2 . (c) I] for g = 0.425 (upper panel) and g = 0.426 (lower 
panel, period doubled LC spectrum). For visualization, a 
small artificial Lorentzian line width was added to all emission 
lines. 

In the fixed points F ± n\^ = <)> = 0, and the time 
dependence comes from the evolution of the total phase 

= —Clot, 

Mt) = ( 8 ) 

with constants C^. Thus the condensates emit light at 
frequency CIq fully synchronized. The emission spectrum 
consists of only one peak, in contrast to the case of non¬ 
interacting polaritons, where two separated peaks are ex¬ 
pected. 

Since the limit cycles LC 1 * 1 are characterized by an 
equidistant spectrum, we first numerically compute the 
corresponding frequency positions, and then calculate the 
intensity of each frequency harmonics using a Fourier 
series expansion. The resulting spectra are shown in 
Fig(2ja,b) . Close to the Hopf bifurcation, there is only 
one considerable emission peak originating from the F + 
spectral line (Fig{2](a,b) upper panels). Further away 
from the Hopf bifurcation, the satellite peaks grow to 
form a frequency comb with asymmetric tails (Fi gHa,b) 
lower panels). The comb also acquires several peak max¬ 


ima, with the highest peak originating from a satellite 
with nonzero N = 2 (Fig{2](a,b) lower panels). When the 
LC undergoes a period doubling bifurcation, the comb 
becomes twice as dense (Figj2](c)). 

The typical modulation frequency is independent of 
the polariton-polariton interaction constant a and is 
of the order of the coupling constant A Cl ~ 7 . For 
typical dissipative rates in semiconductor microcavities 
r ~ 10 12 ps -1 one expects condensate pairs with dissi¬ 
pative coupling 7 ~ 5 x 10 11 ps -1 and below. This would 
permit to generate frequency combs with terahertz sep¬ 
aration between the individual peaks. Additional reduc¬ 
tion of this separation by period doubling can shift the 
modulation frequency into the millimeter range. 

Finally, we consider the influence of noise in Eq. § 
In general, it will broaden the peaks discussed so far, 
and can lead to a merging of peaks with small enough 
spacing. The emission spectrum can be obtained using 
the Wiener-Khinchin theorem, 

OO 

4(H) = hstj (9) 

0 

where ( 0 )) is the auto-correlation of the, now, 

random process i/>n(t). 

The F ± states are periodic orbits in the full four¬ 
dimensional phase space, and the dynamics along these 
periodic orbits is parameterized by the total phase < f>. 
While fluctuations off the periodic orbit will relax back, 
fluctuations along the orbit do not, and will enforce dif¬ 
fusion of $ on the orbit. The latter fluctuations can 
be shown to form a Lorentzian line with the full width 
at half maximum (FWHM) given by W(ni + ri 2 )/ 8 nin 2 
[22] . Note that the FWHM is inversely proportional to 
the number of particles in the condensate, as it should 
be for a laser. 

In contrast to the F ± states, the LC 1 * 1 states are formed 
by the motion on a torus in the full phase space. The 
stability of the attractor demands that fluctuations off 
the torus relax back. Fluctuations along the torus surface 
enforce a diffusion on it. The two nontrivial phases which 
diffuse, are the total phase $ and the second phase angle, 
which characterizes the position on the limit cycle. 

Close to the Hopf bifurcation, and in the presence of 
only a few satellite peaks, we can obtain a closed formula 
for the line width. To parameterize the LC we introduce 
two time arguments: one originating from the total phase 
and the other from the LC phase. Noise in these time 
arguments, according to Eq. ([ 6 ]), leads to 

Mt) =9p(t+ J P(r)dr^j t+ifF^r 

— ^ C N e~ iNA ( t+ ^C) f p ( T ) dT ) e - in ° t e i f F ( T ) dT 1 

N 


(10) 
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where the periodic function p(t) has been expanded in 
a Fourier series with coefficients C N and v(t) is the ve¬ 
locity of the noise-free trajectory along the LC in the 
three-dimensional space {ni,n, 2 , <p}- The noise term 
P(t) is the projected noise along the LC, while F(t) = 


(l/4i) 


- c -c. 


is the noise added to the 


rhs of Eq. ^ for 4>. Note that perturbing the time ar¬ 
gument of p(t) already accounts for a part of the noise 
F(t), because $(t) is time periodic. However close to the 
Hopf bifurcation this periodic part is negligible compared 
to the DC part of 4>(f), so that the above separation of 
noise in the time arguments is justified. 

The choice of the two time arguments is convenient 
when characterizing the dynamics of the noise-free LC. 
In the presence of noise however it leads to nonzero cor¬ 
relations 


((NP(t)+F(t))(NP(t')+F(t'))) = 2(11) 

Denoting ( P(t)P(t ')) = KppS(t — t'), (F(t)F(t')) = 
KFF5(t — t'), ( P(t)F(t ')) = KpF<)(t — t'), we have 

kjv = kff T 2Nkpf + N~npp (12) 

with an asymmetric dependence of km on N. 

The intensity of the noise /cjv(i) is periodic in time, 
which originates from the oscillation of occupation num¬ 
bers and relative phase for evolution along the LC. Ex¬ 
perimentally, the measurement time spans many LC 
periods and one can use the average value Rn = 

T 

4t / kn(t)cIt. Then 
o 


««) = Fi £ l c . 


iN\ 

> I 


kn 


N 


R 2 N + (n-n 0 - NAri) 


(13) 


We obtain a Lorentzian for every emission line, with 
an IV-dependent width Rn according to Eq. (11). The 
^-dependence of the line broadening Rn follows from 
Eq.(12) and shows two remarkable features. First, there 
is an IV-symmetric line broadening with increasing N 
which is very strong and proportional to N 2 . Second, 
there is an asymmetric contribution ~ N which origi¬ 
nates form nonzero correlations kpf- It may lead to a 
satellite peak becoming more narrow than the main peak 
IV = 0, and can further enhance the asymmetry of the 
spectrum, as compared to the noise-free case. 

To calculate the line width R numerically, we denote by 
(A(t), B(t),C(t)) T the normalized tangent vector along 
the LC in the coordinates (0, n\, « 2 ) T - Then P is the 
noise in </>, hi, fi 2 projected onto this tangent and we can 
evaluate from Eq. © 


Wi 


k n = — 

Til 




+ B 2 n 



W± 

n 2 





(14) 
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FIG. 3: (Color online) Noisy (green lines) against noise-free 
(black frequency combs) spectrum of a LC + at g = 0.7, a; = 
0, J = 0.5, for (a) condensation center 1 and (b) condensation 
center 2. The arrows indicate the central peak. Here W = 
0 . 02 . 


We show the spectrum with inhomogeneous line broaden¬ 
ing compared to the noise-free case in Fig(3] Due to noise, 
the asymmetry of the spectrum is enhanced and the strict 
equidistance of emission lines is relaxed for strong enough 
line broadening. We note that this shape together with 
the equidistance of emission peaks is very reminiscent of 
experimentally obtained spectra in [7]. 

Dissipative coupling between coexisting exciton- 
polariton condensates in semiconductor microcavities to¬ 
gether with strong polariton-polariton repulsion leads to 
a rich dissipative nonlinear dynamics already for two cou¬ 
pled condensates. We showed that, in addition to full 
synchronization mm, formation of limit cycles gives rise 
to frequency combs of equidistant asymmetric spectral 
lines. The frequency offset and line spacing of the combs 
are tunable through the control parameters. Through 
period doubling, the line spacing can be additionally re¬ 
duced by an order of magnitude. This modulated emis¬ 
sion can be useful for terahertz and sub-terahertz appli¬ 
cations. Shot noise from the pump results in a complex 
diffusion in phase space and has strong impact on higher 
order satellite peaks. 
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SUPPLEMENTAL MATERIAL 


$ = —Hq- The presence of noise gives 


This Supplemental Material presents the derivation of 
the noise-induced line properties of the F ± limit cycles. 
Let us rewrite the order parameters from 0 as 

Vd ,2 = v/ru^ 1 ’ 2 . (SI) 


4> = -n„(i) + F(t) , F(t) = + F 2 (t )) . (S5) 

Since hi t2 = = 0 we obtain 


From the equations of motion <[2j) it is straightforward to 
find the functions A, B of the differential equations (we 
omit these details for brevity) 


(F(t)F(t')) = 2n5(t — t'), k 


W 

16 


1 

ni 


1 

n 2 


(S6) 


ni ,2 = Ay 2 ( 711 , 2 , <£> 1 , 2 ) , <£ 1,2 = -By 2 (^ 1 , 2 , ^ 1 , 2 ) • (S2) The emission line is finally given by 


Adding noise back to the equations of motion 0 yields 

^ = B,, 2 ( t ) + F M ( t ),F„=i(A- |). (S3) 

The random process F ;j (f) is characterized by 

IT 

(F f j.(t)F ll (t )) = 2 Kn{t)5{t - t’) , p . (S4) 

In a limit cycle the noise-free evolution of tpi )2 {t) is 
determined by the evolution of the total phase given by 




-R / di<e ° 


if F(r)dr 




)4 > 


J df(e- K V (n -° o)t ) 


(S7) 


0 

_ k _ 

7T K 2 + (fl — fig) 2 ! 


which is a Lorentzian centred at f lg and of width k. 






